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Abstract 

The forward-backward asymmetries of the processes e + e~ — > Z — > bb and 
e + e~ — > Z — > cc were measured from a sample of hadronic Z decays col- 
lected by the DELPHI experiment between 1993 and 1995. Enriched samples 
of bb and cc events were obtained using lifetime information. The tagging of 
b and c quarks in these samples was based on the semileptonic decay channels 
b/c — > X + /i and 6/c — > X + e combined with charge flow information from 
the hemisphere opposite to the lepton. 

Combining the and A^ B measurements presented in this paper with pub- 
lished results based on 1991 and 1992 DELPHI data samples, the following pole 
asymmetries were obtained: 

= 0.1021 ± 0.0052 (stat) ± 0.0024 (syst) 
A% = 0.0728 ± 0.0086 (stat) ± 0.0063 (syst) 

The effective value of the weak mixing angle derived from these measurements is 

sin 2 fl^Sff = 0.23170 ± 0.00097. 
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1 Introduction 

The polar angle, 9, of the final state fermion relative to the incoming electron in the 
reaction e + e~ — > // , at y/s ~ M z , is distributed according to: 

dc 8 

oc 1 + cos 2 6> + — Ap B cos6> . (1) 



dcos# 3 

The coefficient of the cos6 term, in the Electroweak Standard Model and for pure 
Z exchange, is related, at the lowest order, to the vector (vf) and axial vector (a/) 
couplings of the Z to the fermions by: 



Afi 3 , a _ 3 2a eVe 2a f Vf 

4 4 a* + t;2 a 2 , + u J 



Higher order electroweak corrections can be accounted for in the above relations by defin- 
ing the modified couplings Vf and df and an effective value sin 2 0^ cff of the weak mixing 
angle: 

^ = l-4| g/ |sin 2 ^ efr (3) 
a f 

where qj is the electric charge of the fermion in units of the proton charge. The effective 
value of the weak mixing angle estimated in this paper is the one corresponding to the 
leptons (sin 2 # w p * ff ), small contributions specific to the quark sector being corrected for 
using the program ZFITTER [1]. 

Because of the values of the Z couplings to fermions, both the forward-backward 
asymmetry and its sensitivity to sin 2 Oyj'lfi are l &r g er i n the Z — > qq channel than in the 
leptonic ones, thus making the Ap^ and Ap B measurements of particular interest. The 
determination of the quark asymmetries Ap B and Ap B requires: 

- the tagging of the Z boson hadronic decays into bb and cc heavy quark final states; 

- the reconstruction of the polar angle of the produced quark/ ant i- quark axis; 

- the orientation of the corresponding axis as a function of the quark direction 1 . 

The analysis presented here is based on events with identified muons or electrons 
produced in semileptonic decays of b and c hadrons, referred to as the "lepton sample" 
in the following. The main parameters used to analyse these events are: 

- the kinematic variables associated with the lepton, namely the transverse {jpr) and 
longitudinal {p£) momentum with respect to the direction of the closest jet; 

- the sign of the lepton electric charge. 

Prompt leptons with high pt and pi allow the selection of a high purity sample of 
e + e _ — > Z — > bb events and, at the same time, the discrimination between quark and 
anti-quark jets on the basis of the charge correlation between the lepton and the parent 
quark. Decay chains like b — > c — > Z + and B°B° mixing reduce this charge correlation. 
Conversely the presence of background and the reduced charge correlations limit the use 
of the largest fraction of the lepton sample at low p T and p L . Two additional variables 
were used in the present analysis to overcome these limitations in the Ap^ measurement: 

- a 6-tagging variable, based mainly on the probability to observe a given event, 
assuming the tracks come from the primary vertex, to isolate pure samples of 
e + e~ — > Z — > bb events; 



J This requirement implies that jets induced by a quark or by an anti-quark have to be distinguished. 
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- a momentum weighted average of the particle charges in the hemisphere opposite to 
the lepton, to provide an independent estimator of the charge of the primary quark. 

By combining the information from the 6-tagging and the lepton pi and pt, a clean 
sample of Z — > cc could also be selected, allowing the measurement of A^ B . 

The thrust axis (T) of the event, oriented by the jet containing the lepton, was used 
to determine the direction of the primary quark. 

The choice of these variables was driven not only by the objective of optimising the 
statistical precision of the measured asymmetries but also by the capability of calibrating 
them on the data, thus controlling well the systematics. 

The data used here were collected between 1993 and 1995 at energies around the 
Z peak with the DELPHI detector at LEP. This analysis extends the previously published 
results based on the events collected in 1990 [2], 1991 and 1992 [3]. 

After a brief presentation of the DELPHI detector, the event and lepton selections 
are described. The observables used in the analysis are discussed together with their 
description by the simulation and the associated sources of systematics. The measurement 
of the asymmetries and A^ B is presented in the last sections. 



2 Detector description and event selection 

2.1 The DELPHI detector 

The DELPHI detector has been described in detail in [4]. Only the components which 
were most relevant to the present analysis are discussed. 

The innermost detector in DELPHI was the Vertex Detector (VD), located just outside 
the LEP beam pipe. It consisted of three concentric cylindrical layers of silicon microstrip 
detectors at average radii of 6.3, 9.0 and 10.9 cm from the beam line, called the Closer, 
Inner and Outer layer, respectively. During 1993 it provided only the measurement of the 
_R$ 2 coordinate and the polar angle acceptance for a particle crossing all the three layers 
was limited by the extent of the Outer layer to 44° < 9 < 136° [5]. In 1994 the Closer and 
the Outer layers were equipped with double sided silicon detectors, also measuring the z 
coordinate [6] . At the same time the angular acceptance of the Closer layer was enlarged 
from 30° < 9 < 150° to 25° < 9 < 155°. The measured intrinsic precision was about 
8 fim for the _R$ measurement while for z it depended on the polar angle of the incident 
track, going from about 10 /im for tracks perpendicular to the modules, to 20 /iin for 
tracks with a polar angle of 25°. For charged particle tracks with hits in all three _R$ 
VD layers, the impact parameter 3 precision was ctr$ = [61/(psin 3//2 9) © 20]yum while for 
tracks with hits in both the Rz layers it was a z = [67/ (psin 5 / 2 9) © 33]/im, where p is the 
momentum in GeV/c . 

Outside the VD, between radii of 12 cm and 28 cm, the Inner Detector (ID) was 
located, which included a jet chamber providing up to 24 /?$ measurements and five layers 
of proportional chambers with both i?<3> and z information. The ID covered the angular 
range 29° < 9 < 151°. In 1995 a new ID was operational, with the same wire configuration 
in the inner drift chamber but a wider polar angle acceptance of 15° < 9 < 165°. 

The VD and the ID were surrounded by the main DELPHI tracking device, the Time 
Projection Chamber (TPC), a cylinder of 3 m length, of 30 cm inner radius and of 122 cm 

2 In the DELPHI coordinate system, z is along the direction of the incoming electron beam, 3> and R are the azimuthal 
angle and radius in the xy plane, and 8 is the polar angle with respect to the z axis. 

3 The R<I> impact parameter is defined as the distance between the point of closest approach of a charged particle in the 
xy plane to the reconstructed primary vertex. The distance in z between this point on the charged particle trajectory and 
the primary vertex is called the z impact parameter. 
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outer radius. The ionisation charge produced by particles crossing the TPC volume was 
drifted to the ends of the detector where it was measured in a proportional counter. Up 
to 16 space points could be measured in the angular region 39° < 9 < 141°. The analysis 
of the pulse height of the signals of up to 192 sense wires of the proportional chambers 
allowed the determination for charged particles of the specific energy loss, dE/dX, which 
was used for particle identification. 

The Outer Detector (OD) was located between radii of 198 cm and 206 cm and con- 
sisted of five layers of drift cells. 

In the forward regions two sets of planar wire chambers, at ± 160 cm and ± 270 cm 
in z, completed the charged particle reconstruction at low angle. 

The muon identification relied mainly on the muon chambers, a set of drift chambers 
providing three-dimensional information situated at the periphery of DELPHI after ap- 
proximately 1 m of iron. One set of chambers was located 20 cm before the end of the 
hadronic calorimeter (HCAL), two further sets of chambers being outside. At 9 ~ 90° 
there were 7.5 absorption lengths between the interaction point and the last muon detec- 
tor. 

In the Barrel part of the detector there were three layers each including two active 
planes of chambers covering the region | cos#| < 0.63. The two external layers overlapped 
in azimuth to avoid dead spaces. In the Forward part, the inner and the outer layers 
consisted of two planes of drift chambers with anode wires crossed at right angles. The 
resolution was 1.0 cm in z and 0.2 cm (0.4 cm) in for the Barrel (Forward). In 1994 
a further set of chambers (Surround Muon Chambers) was added to cover the region 
between the Barrel and Forward chambers. 

The electromagnetic calorimeter in the barrel region, |cos#| < 0.73, was the High 
density Projection Chamber (HPC), situated inside the superconducting coil. The detec- 
tor had a thickness of 17.5 radiation lengths and consisted of 144 modules arranged in 6 
rings along z, each module was divided into 9 drift layers separated by lead. It provided 
three-dimensional shower reconstruction. In the forward region, 0.80 < |cos#| < 0.98, 
the electromagnetic calorimeter EMF consisted of two disks of 5 m diameter with a total 
of 9064 lead-glass blocks in the form of truncated pyramids, arranged almost to point 
towards the interaction region. 

2.2 Selection of hadronic events 

The selection of charged particle tracks and neutral clusters was performed according 
to the requirements of Table 1. Hadronic events were then selected with an efficiency of 
95% requiring: 

• at least 7 accepted charged particles; 

• a total measured energy of these charged particles, assuming pion masses, larger 
that 0.15 times the centre-of-mass energy, s/s. 

A total of 2.7 million hadronic events was selected from 1993-95 data, at centre-of- 
mass energies within ± 2 GeV of the Z resonance mass. A set of 8.4 million simulated 
hadronic events for years 1993 to 1995 was used, generated using the JETSET 7.3 Parton 
Shower model [7] in combination with the full simulation of the DELPHI detector. The 
parameters of the generator were tuned to the DELPHI data as described in [8]. The 
detailed breakdown of the events used in data and simulation for each year is given in 
Table 2. 
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charged-particle tracks 


polar angle cos 9\ 


< 0.93 




length of track measured inside TPC 


> 30 cm 




impact parameter (R&) 


< 5 cm 




impact parameter \z\ 


< 10 cm 




charged particle momentum 


> 0.2 GeV/c 




relative uncertainty on the momentum 


< 100% 


neutral clusters 


detected by HPC or EMF 






polar angle cos# 


< 0.98 




HPC (EMF) energy 


> 0.8(0.4) GeV 



Table 1: requirements on charged particle tracks and neutral clusters for hadronic events 
selection. On selected events, the energy flow measurement and jets reconstruction 
were performed with an improved neutral clusters reconstruction, including photons of 
lower energy in the HPC (down to ~ 0.3 GeV ) and neutral showers of more than 1 
GeV reconstructed in the HCAL. 



Year 


# of events (in 10 3 ) 




Data 


Simulation 


1993 


696 


2276 


1994 


1370 


4300 


1995 


662 


1829 



Table 2: The number of selected hadronic events for data and simulation 

3 Lepton samples 

The main kinematical variable used to measure the flavour composition of the leptonic 
sample was the transverse momentum, p T , of the lepton with respect to the jet axis. To 
compute this axis the jet containing the lepton was used, but its direction was recon- 
structed without the lepton. Jets were built using the JADE algorithm [9] with a scaled 

2 

invariant mass y cut = -^r~ > 0.01 . 

vis 

To ensure a good determination of the thrust polar angle, $t, the analysis was limited 
to events with looser)! < 0.9. For events with more than one lepton candidate, only the 
highest pt lepton was used for the Apg and A^ B measurement. 

The lepton identification has been studied not only by means of special data and 
simulation samples (for example : K° — > 7r + 7r~ , Compton events) but also using 

p, pt and 6-tagging 4 cuts to select, in hadronic events, lepton sub-samples with different 
purity. In practice with only two such lepton sub-samples, two parameters, the efficiency 
and the purity of the overall sample, can be compared between data and simulation. 

4 The 6-tagging method and its calibration will be presented in detail in Section 4. 6-tagged and anti-6-tagged samples 
refer, respectively, to the purer and to the most contaminated 6 sub-samples as defined in Section 4. The 6-tagging is used 
in this section to estimate the purity of the electron sample but is only used as a cross-check for the muon sample. For this 
reason to avoid any sizable correlation in the tuning of the simulation, the 6-tagging calibration used in this section was 
based on events from the muon sample alone. 
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Year 


Energy 


muons 


electrons 




89 GeV 


6068 


4240 


1993 


91 GeV 


28791 


21553 




93 GeV 


9171 


6536 


1994 


91 GeV 


95183 


69971 




89 GeV 


5147 


4062 


1995 


91 GeV 


28600 


21786 




93 GeV 


8629 


6443 



Table 3: Number of events with at least one lepton candidate for the different years and 
for the three centre- of-mass energies. The highest Pt lepton is used to classify the event 
as "muon" or "electron". 

The number of lepton candidates for the different years and centre-of-mass energies can 
be found in Table 3. Details on the lepton identification and on the sample composition 
are given in the next two sub-sections. 

3.1 Muon sample 

For the muon identification the tracks reconstructed in the central detectors were 
used to define a path along which hits in the muon chambers were looked for. The 
identification algorithm has been described extensively in ref . [3] . Muon candidates with 
a momentum, p, above 2.5 GeV/c and in the region of good geometrical acceptance were 
selected. The muon polar angle 9^ was required to be in the region 0.03 < |cos6>^| < 0.6 
or 0.68 < | cos ^ | < 0.93. Only for a small fraction of the 1994 data sample the Surround 
Muon Chambers, which filled the gap between the barrel and forward detectors, were 
able to provide useful muon identification. 

The muon identification efficiency was measured in Z — > fi + fi~, Z — > t + t~ and 
77 — » fi + fi~ events, yielding on average about 0.85 for 45 GeV/c muons and 0.83 for 
momenta between 10 and 5 GeV/c . 

In order to extract Apg and Ap B from the observed asymmetry the absolute lepton 
efficiency is not required, only a correct description of the sample purity is needed. The 
contamination from misidentified hadrons arises partly from the decay of pions and kaons 
and, for momenta above 3 GeV/c, mostly from energetic hadrons interacting at the end 
of the calorimeter and generating punch-throughs. Kg particles decaying into two pions 
were used to measure the rate of pion misidentification above 3 GeV/c showing that 
the fraction of pions misidentified as muons was 1.79 ± 0.09 ± 0.05 and 1.41 ± 0.10 ± 
0.03 times bigger in data than in the simulation for the barrel and the forward regions 
respectively. Muons from pion decays were subtracted from the misidentified sample to 
compute the numbers quoted above. The first error is due to the limited statistics, the 
second corresponds to a 15% change in the contamination of muons from pion decays. 
The fraction of muons from tt and K decays present in the selected muon sample was 
determined to this precision by comparing the size of two muon sub-samples in hadronic 
events selected by p,pt cuts 5 . 

To measure further the sample composition directly from the data, the number of muon 
candidates, normalised to the number of hadronic Z decays, was compared between data 
and simulation in sub-samples enriched in prompt muons or background by different sets 

5 Thc first one was selected with p > 2.5 GeV/c and pj> < 0.7 GeV/c and the second one with p > 4 GeV/c and pt > 0.7 
GeV/c , the two samples having respectively 19% and 4% of muons from it and K decays. 
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Lepton candidate source 


Sample com 


Dosition in % 

e 


Lepton from b hadron decay : "same sign" 


32.5 


35.1 


a) b — > l~ 

b) b^T^l- 

c) b — ► c — > l~ ] b — > c — > t — > l~ 


29.0 
1.0 
2.5 


31.6 
1.0 
2.5 


Lepton from b hadron decay : "opposite sign" 


11.8 


11.7 


d) b -> c -> 1+ ; b -> c -> r -> /+ 






Lepton from 6 hadron decay : other source 


0.3 


0.3 


e) b -> J/* -> /+/- 






Lepton from c decay 


16.8 


15.9 


f) c -> /+ ; c -> r -> 1+ 






Background 


38.6 


37.0 


g) Misidentification 

h) Light mesons decay / converted gammas / other 


26.3 
12.3 


18.2 
18.8 



Table 4: Full 1993-1995 lepton sample composition. The leptons from heavy flavour 
decays, when the heavy flavour quark was coming from a gluon splitting, are counted in 
"other" , line "h)". The total efficiency to select a muon or an electron from the process 
"a) "is respectively (44.7 ± 0.2)% and (35.4 ±0.4,)% including all effects (momentum cuts 
and detector inefficiencies) except the efficiency to select hadronic events. 



of selections in p, px or 6-tagging. The stability of the misidentification correction quoted 
above was measured as a function of p , pr and 6 in b and anti-6-tagged samples. No 
discrepancy was observed outside expected statistical fluctuations. For example in the 
anti-6-tagged sub-sample, which had a purity of 30% in leptons from heavy flavour decay, 
a data/simulation comparison for 28 bins in momentum had a 27% x 2 probability. In 
this same sub-sample the 6.6% excess of positive particles due to the difference between 
the K + and K~ cross-sections in the detector, was also perfectly described with a statis- 
tical precision of ± 0.5%. These studies confirmed the correct description in the tuned 
simulation of the different sources of background within the uncertainties quoted above. 

The comparison between data and simulation for the cos 9^ distribution is presented 
in Figure 1. The muon sample composition is given in Table 4. 

3.2 Electron sample 

The electron candidates, of momentum higher than 2 GeV/c, were identified in the 
barrel (0.03 < |cos# e | < 0.7 ) by combining the electromagnetic shower information from 
the HPC and the track ionisation measured by the TPC, with a neural network. In 
the forward region (0.7 < |cos# e | < 0.9) only the ionisation measured by the TPC was 
used. Mainly due to the large amount of material in front of the EMF, the calorimetric 
information in the forward regions was not used to identify electrons. 

Electrons in the barrel region 

In the barrel the contamination and efficiency of the electron sample was tuned in the 
simulation using two sub-samples with p > 3 GeV/c and a b or anti-6-tag. Their content 



7 



of misidentified electrons was respectively 27% and 89%. The efficiency was found to 
be correctly described by the simulation, in agreement with a study based on Compton 
scattering and photon conversion samples [10]. The misidentification probability was 
found in the data to be a factor 0.9 lower than that in the simulation, with variations 
within a few percent as a function of the year and angular region. A study based on 
sub-samples selected by p and pt cuts gave compatible results. The relative precision on 
this correction was estimated to be ±5%. 

The number of selected electron candidates in a low p (p < 3 GeV/c ) and px (j>t < 
1 GeV/c ) region for anti-6-tagged events agreed between data and simulation for the 
standard and low converted gamma rejections (with converted gamma rejections of 90% 
and 75% respectively). This study was statistically compatible with a correct description 
of the converted gamma content of the electron sample at a ± 10% level. 

Electrons in the Forward regions 

In the forward regions, where only the ionisation measured by the TPC was used, 
the misidentification could be studied with the muon sample, muons and pions having 
almost the same ionisation signature in the TPC. This showed the need to increase the 
misidentification by a factor 1.12 ± 0.07 in the simulation. The amount of electron from 
converted photons was determined with the help of the 6-tagging and p, p T requirements 
as above. This amount was found to be correctly described by the simulation for the 
different years within a precision better than ± 10%. 

The comparison between the data and the simulation for the cos# e distribution is 
presented in Figure 1. The composition of the electron sample is quoted in Table 4. 

4 Use of the 6-tagging information 

To improve the separation between heavy and light flavours a 6-tagging technique 
developed for the measurement of Rb, the partial width of Z into bb pairs, was used [11]. 
Each event was divided into two hemispheres according to the direction of the thrust axis, 
and the 6-tagging probability for a hemisphere to contain a b quark was given by the jet 
with the highest 6-tagging probability in the hemisphere. 

The tagging technique was based on the combination of up to four discriminating 
variables, %i (i=l to 4), defined for each jet: 

• the jet lifetime probability, constructed from the positively signed impact parameters 
of all tracks included in a jet, is built from the probabilities to observe a given set 
of impact parameters, assuming the tracks come from the primary vertex; 

• the effective mass of the system of particles assigned to the secondary vertex 6 ; 

• the rapidity of tracks associated to the secondary vertex with respect to the jet 
direction; 

• the fraction of the jet energy carried by charged particles from the secondary vertex. 

The correct description of these variables by the simulation is shown in [11]. 

For each jet containing a secondary vertex, the four variables Xi were then combined 
into a single tagging variable y by means of their probability density functions ff ds {xi), 

6 A secondary vertex was required to contain at least 2 tracks not compatible with the primary vertex and to have 
L I ctl > 4 where L is the distance from the primary to the secondary vertex and ox is its error. Each track assigned to 
the secondary vertex should have at least one measurement in the VD and at least 2 tracks should have measurements in 
both R — <E> and R — z planes of the VD. 
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Figure 1: cos 9 distribution for muon, upper plots, and electron, lower plots, candidates. 
The lepton candidates have px < 0.7 GeV/c in a),c) and px > 0.7 GeV/c in b),d) corre- 
sponding to samples enriched in background and signal respectively. The "b — > n/e right 
(wrong)" samples correspond to leptons with the same (opposite) sign as the primary b 
parton (see Section 6.1). Converted photons cause the strong cos9 e dependence of the 
background in the c) plot. Only the statistical errors are quoted in these plots. The sys- 
tematics on the samples composition considered in the text are enough to cover the few 
discrepancies observed between the data and the simulation. 
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fi(xi), fi(xi), for uds, c and b quarks respectively, as determined from simulation studies: 



y = n, 



n$s+»-ir 



n u ds being the fractions of c-jets and uds-jets with a reconstructed secondary vertex 
normalised by the relation: n c + n u d s = 1. Only the first variable was used if no secondary 
vertex was reconstructed. Hemispheres with a jet containing a b-quark were characterised 
by a large value of the variable t/hem = — log 10 V- 

The value of the tagging variable for the whole event was computed from the corre- 
sponding values obtained in each hemisphere as: 

TjEVT = max(r)HEMl,r)HEM2)- 

The sample composition determination, as a function of the value of 7]evt, was needed 
to extract Ap^ and Ap B from the raw asymmetries. Since a separate tag for each hemi- 
sphere was used, the sample composition could be derived from the data themselves with 
minimum input from the simulation by using a technique similar to the single tag versus 
double tag method of the R b analysis [1 1] . 

For events with the thrust axis situated within the barrel acceptance (|cos#r| < 0.7), 
the distribution of the hemisphere 6-tagging variable t\hem was divided into three intervals 
corresponding, respectively, to events enriched in uds (I), c (II) or b (III) flavours. For 
this low number of intervals a direct measurement of their content in term of uds, c and 
b can be implemented in the data as follows. 

In each interval j the fraction of events with at least one hemisphere in that 
interval and the fraction fjp of hemispheres in the interval itself, were expressed by the 
following relations: 

/if> = ^r q (2e®-e%)=Y,rdn*-P?-%Hl-pW 
<? <? 

ff = w 

q 

where e q ^ were the fractions of hemispheres in the j— th interval for the flavour q (q = 
uds,c,b), and the correlations p q j ^ = (e^ — e q ^ 2 )/(e q \l — e q ^)) accounted for the 

probability ( ) of having both hemispheres in that interval. The variables r q stand 
for the fractions of Z — > qq events in the selected leptonic sample. 

The requirement of an identified lepton in the final state strongly enhanced the fraction 
of events with a Z decaying into heavy quark pairs. Therefore the fractions r q were 
obtained from R q , the Standard Model partial decay widths of the Z , via the relation 

r q = R q — — q = uds,c,b (5) 

Chad/ 

where was the flavour dependent hadronic selection efficiency, taken from the simu- 
lation, and e had/ = J2uds,c,b R q e i/- 

To determine the fractions e^ RD in real data for the different intervals, it has been as- 
sumed that they differ only slightly from the ones in the simulation, e q ^\ 



Jj) _ M) | 

t q,RD — L q 1 



l + #>) 
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In the approximation, confirmed by the data, of small corrections 8 q \ the set of 
Equations (4), including the closure relations on the fractions e^ RD and e q \ gives at first 
order in 8q^: 

J2 r q eP [2 - p& - 2ef (1 - pf )] 8^ = /« - £ r q e^ [2 - pf - ef (1 - p<f > )] (6) 

9 9 

E^4 3) 4 J) = /^-E^4 3) (7) 

<? 9 



where p q ^ and e q ^ were taken from the simulation. For N int = 3 5-tagging intervals there 
are in total 9 unknowns 8 q ^ and 9 equations. The rank of the matrix of the coefficients 
is 8 so that one input 8 q ^ was required 7 . For Ni nt = 2 and merging together 2 flavours, 
the system reduces to 6 equations and 4 unknowns. Since the rank of the matrix of the 
coefficients is 4 the system has one exact solution 8 . Therefore the 5^ for N int = 3 were 
obtained in two steps. First we combined together the two highest bins of r] HEM , merged 
the c and b contributions and solved this reduced system with 4 unknowns. Then the full 
system was solved fixing 8^J S to the value obtained in the previous step. 



For events with 0.7 < |cos 6 T \ < 0.85, because of the reduced performances of 6-tagging 
in the forward region, only the reduced system with 4 unknowns was solved. No 6-tagging 
information was used for events with (cos^t) > 0.85 9 . 

Errors on 8 q ^ due to the finite statistics of the simulated sample were estimated in 
the following way. For each flavour q, we considered the two dimensional distributions 
{Vhemi, VHEM2} which could be derived from the original one in the simulation by adding 
-1, 0, and +1 standard deviations to the content of each interval. This was done conserv- 
ing the total number of events of that flavour and with the standard deviations given by 
the multinomial distribution. For each configuration the coefficients e q \ p q ^ in Equation 
(4) were recomputed and then the system solved. The spread of the different solutions 
for the 5q^ was considered as the simulation statistical error on these corrections. 
As a cross check of the method, the simulated sample was divided into 6 different sub- 
samples of equal size. For each sub-sample, the system was solved and the uncertainty on 
the solutions was evaluated by using the procedure described above. The spread of the 
solutions in the subsets was found to be in agreement with the estimation of the error. 
The corrections 1 + 8 q ^ to the simulation fractions found for the 1994 sample together 
with the error due to the finite simulation statistics are shown in Table 5. For all the 
samples, the corrections S^IP were found to be compatible with zero indicating a good 
control of the background level in the region most relevant in this measurement. The 
fractions e^ 11 ^ were found instead 2-4 % higher in the data than in the simulation. 

For the system with N int = 3 the predicted correlations have a sizable value only for 
p[ IH) ( = 0.027 ± 0.005 in 1994 ). The detector and QCD origins of such correlations 
have been studied in detail in [11]. In the present analysis even a 100% change in the 

7 As by construction X^" 1 * £< q^ = 1 implies X^ int fjj^ = * > onc equation among the Equations (7) and (8) can be 
deduced from the other. 

8 For Ni n t = 2 , we have by definition p q ^ = p^ 2 ' which makes in this case the two equations associated to (6) equivalent. 
9 For the 1993 data sample, due to the reduced length of the micro-vertex detector, the b-tagging was performed only 
down to |cos#t| < 0.81. 
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bin 


i 


II 


III 


(dominant flavour) 




(c) 


(b) 


Barrel region 


M) 

uds 


.71 


.23 


.06 


1 + °uds 


1.019 ± 0.003 


0.950 ±0.017 


0.986 ± 0.072 




.45 


.34 


.21 




0.986 ± 0.007 


1.046 ± 0.013 


0.956 ± 0.031 


Jj) 


.13 


.18 


.69 


1 + 5? 


0.970 ± 0.009 


0.946 ± 0.006 


1.020 ± 0.002 


Forward regions 


U) 

uds 

1 4- 


.70 

1.000 ± 0.009 


.30 

1.000 ± 0.021 


Ji) 

fc fec 


.36 


.64 


i + *2 ) 


.931 ± 0.008 


1.039 ± .005 



Table 5: Values of the fractions ( e q ), obtained from the simulation, and of their 
respective corrections (1 + 5 q ^), fitted on real data, using 1994 event samples. 

predicted correlation has a small impact on the estimated data sample composition. The 
variation induced is of the same order as the one associated to the statistical uncertainty 
on 6®. 

For the system with N int = 2 the predicted correlations were up to ~ 0.1 for the b/c 
flavours and still compatible with zero for the uds sample. The high value of obtained 
in this case is a pure artifact of the merging of the b and c samples and is just related, at 
first order, to the difference in tagging efficiency of b and c. 

The merging of b and c for N int = 2 is justified by the fact that for the 6-tagging 
intervals used in this case, the 8 q ^ corrections are mainly related to the difference in 
the description of the detector response between real data and simulation and not to 
the details of the b and c physics. This is supported by the fact that in the interval / 
dominated by uds, which is the same both for N int = 2 and N int = 3, the corrections 
5^ and 5$p in the barrel are compatible (cf. Table 5). To evaluate possible biases from 
this merging procedure, another system, also with Ni nt = 2, was built starting from the 
original one with N int = 3 but now combining the two lowest bins of tjhem, bins I and 
II, and merging the uds and c flavours. The changes found in and A^ B were taken 
conservatively as systematic errors (cf. Section 7). 

5 Use of the jet charge information. 

The jet charge measured in the event hemisphere, opposite to the lepton, provides 
an additional information on the charge of the parton from which the lepton originates. 
This information is particularly relevant for events with a low pt lepton to be still able 
to distinguish between b — > l~ and b — > c — > / + and for events with a high p T lepton 
to tag B°B° oscillations. The jet charge was built by means of a momentum-weighted 
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(Pi) average of the charges (qi) of the charged particles in the hemisphere opposite to the 
lepton: 

r\ ^hem Qi\Pi ' T\ 

^°PP v ,-> 7ft {K 
l^hem \ Pi ' 1 I 

with the event divided into two hemispheres by a plane perpendicular to the thrust axis 
T. 

With this definition the information coming from the tracks in the lepton hemisphere 
was not used in order to avoid the strong bias in the topology due to the presence of a 
lepton. Based on the work presented in [12] K = 0.8 was chosen to optimise the b/b 
separation. We restricted the use of the jet charge to the events with the thrust axis in 
the barrel ( |cos0t| < 0.7) and belonging to the lepton subsample enriched in b (bin III 
according to the definition given in Section 4). 

The distribution of the total event jet charge in the hadronic decays of the Z turned 
out to be systematically displaced from zero (~ +0.01), due to the hadronic interactions 
of particles inside the detector. It was checked that this shift was independent of the event 
flavour in each 6-tag bin and is corrected for, separately in the data and in the simulation, 
as a function of the thrust axis of the event. After this correction it was possible to treat 
in a consistent way positive and negative leptons, by using, as b/b discriminating variable, 
the product of the lepton charge times the opposite jet charge, Qs> x Q opp . For a pure 
sample of leptons coming from Z — > bb decays, Qe x Q opp has a Gaussian distribution 
centred at negative (positive) values in case of right (wrong) sign correlation between the 
lepton and the b parton from the opposite hemisphere. In the following this central value 
will be quoted as ±Q - et j and the width of the Gaussian as ctq t ( . After normalisation 
by the total number of leptons, the integral of the Gaussian with the negative mean will 
be quoted as /, the integral of the other Gaussian being then equal to 1 — /. 

A procedure of self calibration of Qe x Q opp with the data was used so as not to rely 
on the simulation for the jet charge description. It also allowed the fraction / to be more 
independent of the precise knowledge of the B°B° mixing or of the branching fractions 
for the direct (b — > l~) and cascade (b — > c — > i + ) semileptonic decays. 

The first step of the jet charge self calibration consisted in the tuning of the simulation 
in order to reproduce the total event jet charge distribution measured in data. The total 
event jet charge measured in the full sample of 6-tagged hadronic decays could be used. 
This distribution gives a direct estimate of <jQ jetl [12]. 

As a second step the values of / and Qj et l were obtained by a double Gaussian fit of 
Qe x Qopp in the data. This has been done for each year for muon and electron separately, 
after subtracting the background predicted by the simulation. The statistical sensitivity 
of the fit was improved by reducing the number of fitted parameters using the following 
constraint: 

< Qt x Q opp >= (1 - 2f)Q jetjl . (9) 

This constraint is derived from the definition of the two Gaussian distributions introduced 
above. 

Subsequently the jet charge values in the simulation were corrected to reproduce the mea- 
sured distribution of mean, Qj ett i, and width, ctq . Moreover, in each lepton subsample, 
the events were re-weighted in order to reproduce the fitted value of the fraction of right 
sign leptons, /. After this calibration the simulation describes the data correctly as can 
be seen in Figure 2; event sub-samples enriched, by kinematical cuts, in leptons from 
different origin like b — > l~ in the high p^ region or b — > c — > l + in the low pt one, are 
well described even if they have quite different values for < Q e x Q opp >. Figure 2 gives 



13 




Figure 2: Jet charge distribution for the full muon sample. The muon candidates have 
Pt < 1.3 GeV/c and p < 7 GeV/c (left) and px > 1.3 GeV/c (right) corresponding to 
samples enriched in b — > c — > yU + or 6 — > /i _ events respectively. 

a good consistency check of the overall simulation tuning regarding the lepton sample 
composition. 

There are two sources of uncertainty related to the jet charge self calibration method 
(see Table 6) : 

• the o-Q jet l , / and Qj et) i are measured with a statistical uncertainty in the simulation 
and in the data, 

• the (JQ. etV f and Qj etl values are extracted from the data after subtraction of 
the uds and c contamination. This subtraction, as it relies on estimates from the 
simulation, induces systematic errors. 

These errors have been estimated using samples enriched in uds or c events, corre- 
sponding to the 5-tagged events in bins (I) or (II). In such samples any difference in 
the event charge or jet charge distributions between data and the tuned simulation was 
entirely attributed to an imperfect description in the simulation of the uds or c events. 
Using a new description of the event or jet charge for uds or c events to fix these differ- 
ences, new values of o-Q jetl OY f and Qj et) i were then estimated. The biggest changes in 
a Qjet i ' / an d Qjet,i observed have been considered as systematic errors. 



6 The fit of the asymmetries 

The and A^ B asymmetries were extracted from a minimum x 2 fit t° the observed 



charge asymmetry, A°p^\ defined as: 



obSti _ N-(i)-N+(i) 

FB " ~ 



N-(i) + N+(i) 
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fm% 


Q-ietd 


/i sample 


e sample 


jj, sample 


e sample 


Data 


.2842 ± .0004 ± .0021 


68.4 ± .6 ± .6 


70.2 ± .8 ± .6 


.099 =F -003 =F -003 


.098 =F -004 =f -003 


Sim. 


.2894 ± .0003 


69.3 ± .4 


71.1 ± .4 


.103 =F .002 


.101 =F -002 



Table 6: Values of aQ jetl , f and Qj et i in the barrel for the subsample enriched in b, fitted 
in the data and in the simulation before tuning in 1994- It should be noticed that, because 
of the constraint expressed by Equation (9), f and Qj et j are fully anti- correlated. The 
first quoted error corresponds to the analysed statistics and the other, in the case of data, 
to systematics. 

where N + (i) and N~(i) are the numbers of events with an identified lepton in the i-th 
bin with, respectively, a positive and a negative electrical charge. 

Four variables were used for binning the sample: cos 9t, which accounted for the polar 
angle dependence of the asymmetries, and three multivariate classification parameters, 
chosen to have bins enriched with leptons from a single origin. These last three parameters 
allowed the statistical errors of the and Ap B measurements to be reduced. 

6.1 The multivariate parameters 

The observables entering the multivariate parameters, whose values depend on the 
origin of the lepton candidates, were chosen to be: 

• the transverse (pr) and longitudinal (pi) momenta of the lepton; 

• the event 6-tagging, r] EVT ; 

• the product of the lepton charge times the jet charge of the opposite hemisphere, 

Qi x Qopp- 

Starting from these observables a multivariate tagging of the lepton origin was built 
by considering four classes: 

1. b r : leptons from b hadron decays in Z — > bb events with the right sign correlation 
(same sign) with respect to the primary b quark; 

2. b w : leptons from b hadron decays in Z — > bb events with the wrong sign correlation 
(opposite sign) with respect to the primary b quark; 

3. c : prompt leptons from c decays in Z — > cc; 

4. bkg : other processes (misidentified hadrons, leptons from light hadron decay, elec- 
tron and positron from photon conversion and leptons from heavy flavour hadron 
decays where the heavy flavour quarks were produced by gluon splitting). 

The sign correlation mentioned here refers to the one between the lepton charge and 
the b/b flavour at production and therefore it includes possible effects due to B°B° mixing 
(cf. Section 6.3 for a more extended discussion on the mixing). The probability densities 
p^ T,PL and pt a 9>i et ~ ch f observing a set of (pr, Pl) and (tjevt, Qi x Qopp) values for a lepton 
from the class k were computed by using two-dimensional distributions from the tuned 
simulation. A likelihood ratio Vk was built to estimate the probability corresponding to 
a given set of values within a class: 
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where Nj, (N^r) is the total number of leptons from the class k (k'). The scaling of the 
likelihood ratio by N k takes into account the relative weights of each class. Neglecting 
some of the correlations between the observables, such a definition identifies Vk as the 
fraction of lepton candidates with a given set of Pt,Pl,Vevt and Qe x Q opp belonging to 
the class k. 

This technique, used for the multivariate classification, extends that in [11] by con- 
sidering probabilities in two dimensions and takes into account part of the correlations 
between pairs of observables. 

In order to consider the possible improvement by taking into account all possible corre- 
lations, an approach based on a classification with a neural network was also tried. The 
results obtained were in good agreement with those from the multivariate approach but 
had a slightly worse statistical precision. The multivariate approach was chosen as it 
allowed, for the small number of observables used, a simpler control of the analysis and 
an optimal use of the available simulation statistics. 

The Vb r ,Vb w and V c distributions for muons and electrons are shown in Figure 3. 

6.2 Measurement of A 1 ^ and Af B 

The asymmetries and Ap B were extracted from a x 2 fit to the observed asymme- 
tries A°p S g over the different bins of the (cos#t, Vb r -Vb w , V c ) parameter space: 

2 ^ (( (& - fO A f* + fi A f* + n k A k B l ) w e T - a °fA 2 
x 2 = }^— -2 — ( l0 ) 

i % 

where : 

• Wl = f tth^— Y]^-f a — c ° s % „ takes into account the dependence of the asymmetry 

6t ?> N data^l- X 1 + (C0S4) 2 J J 

on the polar angle; 

• oi is the statistical error including contributions from both data (Ap^ 1 ) and simu- 
lation (f l k ,A b FBl; 

• fl are the different fractions in each bin determined from the tuned simulation. 

To optimise the use of the available statistics, the multivariate variables were computed 
separately for the different years and lepton samples and all samples were merged for the 
X 2 fit. The data binning in the parameter space was done to have the same number of 
events per bin, ~ 100 , ~ 180 and ~ 150 events for = 89.43 GeV , 91.22 GeV and 
92.99 GeV respectively. 

Due to the opposite sign in the contribution of the b r and b w classes to the b asymmetry 
only the difference between the fractions of b r and b w classes matter in practice. It's 
why the equi-populated bins for the x 2 fit were defined by using the combined variable 
Vb r — Vb w and V c . A possible third sampling corresponding to Vbkg has been found of 
marginal interest, mainly due to the closure relation on the Vk, and has not been used 
for the present fit. 

A sign correlation between the lepton candidate and the parent quark can exist also for 
the misidentified leptons thus leading to non-zero values for the background asymmetry 
A b pQ l . Furthermore, since this correlation increases with the particle momentum and 
as a function of 6-tagging value, A b p^ 1 must be known in each bin. To optimise the 
statistical precision of the estimated A b p^\ the same factorisation technique as in the 
previous analysis [3] was adopted: the simulation was only used to determine the charge 
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Figure 3: Likelihood ratio Vb r (up) , Vb w (middle) and V c (bottom) distributions for 
muons (left) and electrons (right) for all years. 
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correlation between the background and the initial quark in each bin, while the quark 
asymmetries were set to their Standard Model expectation for background in light quark 
events 10 or to the fitted parameters , A^ B for background in b or c events. 

6.3 Effect of the B°B° mixing 

The B°B° mixing reduces the charge correlation between the initial b/b produced from 
the Z decay and the lepton issued from the B hadron semileptonic decay. The size of the 
change depends on the proper decay time of the B hadron and on its type, resulting in 
different values of the effective mixing in the different bins of the lepton sample, for the 
following reasons: 

• the B® and B° s fractions in the b sample are not the same for direct or cascade 
decay leptons due to differences between the D + /D° /D s production rates from the 
different B hadrons; this introduces, for example, a variation of the effective mixing 
as a function of Pt] 

• the use of the 6-tagging biases the content of the bins in terms of proper decay time, 
thus introducing sizable changes in the effective mixing; 

• the sign correlation between the lepton and the jet charge, measured in the opposite 
hemisphere, depends directly on the mixing. 

The B°B° mixing is now well measured [13]. Following the approach developed in the 
LEP oscillation working group [14] , the simulation was tuned to reproduce the measured 
B fractions ( f B ± , f B o, f B o ,f'B baryon ) and the time dependence of the oscillations ( Arrid 
and Am s ). The values and the corresponding uncertainties used to implement the B°B° 
mixing in the simulation are listed in Table 7. 11 

With this approach, the values estimated from the tuned simulation of and 
included the expected amount of mixing. 

6.4 Results 

The measured asymmetries and the corresponding statistical errors using the 1993- 
1995 lepton samples are listed below: 

At y/s = 89.43 GeV : 

Af B = 0.066 ± 0.022(stat) 

Af B = 0.030 ± 0.035(stat) 
with a correlation of 0.19 and ^ = ||| , Prob(x 2 )=0.87; 

at ^ = 91.22 GeV : 

A% = 0.0958 ± 0.0061(stat) 

1(1 The uncertainty due to the exact knowledge of these asymmetries is negligible compared with the error on the charge 
correlation itself. 

11 The values used come from the LEP Lifetime Working group (lifetimes), the LEP oscillation working group (fractions, 
Amj and Am s ) and the LEP Heavy Flavour working group (x). The lower bound value quoted for Am s was used 
(Am 3 >10.6 ps _1 ), no sensitivity to the exact value of this parameter in the allowed domain has been observed. All these 
numbers are taken from [13]. 
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Af B = 0.0585 ± 0.0098(stat) 
with a correlation of 0.22 and ^ = ff|§ , Prob(x 2 )=0.75; 

at ^ = 92.99 GeV : 

Af B = 0.109 ± OMS(stat) 

Af B = 0.108 ±0.028(stat) 
with a correlation of 0.19 and ^ = , Prob(x 2 )=0.57. 

The result of the fit for a data subsample is presented with the observed asymmetry 
in Figure 4. 
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Figure 4: Observed A°p s B asymmetry at peak energy as a function of cos9 T for a data 
subsample enriched in b with the right sign (Vb r — Vb w > 0.6). The result of the fit is 
shown as a dashed line. 



7 Systematic effects 

The systematics from the different sources are listed in Table 7. 
Lepton sample 

To estimate the systematics due to uncertainties in decay branching ratios and spectra, 
the standard prescription of the LEP Heavy Flavour Working group was used [15]. The 
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Parameters 


Central 


Variations 


A ^FB 






value 


applied 


Peak 


Peak 


Br(b -> /) 


0.1056 


± 0.0026 


=F 0.00048 


± 0.00062 


Br{b -> c -> Z) 


0.0807 


± 0.0034 


=F 0.00015 


=F 0.00080 


£r(6 -> c -> Z) 


0.0162 


+0.0044 
-0.0036 


± 0.00017 


± 0.00180 


5r(6 -> r -> Z) 


0.00419 


± 0.00055 


=F 0.00001 


± 0.00027 


5r(6 -> J/^ -> Z) 


0.00072 


± 0.00006 


± 0.00005 


± 0.00002 


5r(c -> Z) 


0.0990 


± 0.0037 


± 0.00036 


=F 0.00182 




0.21644 


± 0.00075 


=F 0.00007 


± 0.00007 


^ cc/^ had 


0.1671 


± 0.0048 


± 0.00034 


=F 0.00130 


g^bb 


0.00254 


± 0.00051 


± 0.00012 


± 0.00005 


g -> cc 


0.0296 


± 0.0038 


± 0.00012 


± 0.00001 




0.702 


± 0.008 


=F 0.00016 


=F 0.00019 


(Xe) d * in cc events 


0.510 


± 0.0094 


± 0.00047 


=F 0.00046 


6 decay model [15] 
c decay model [15] 


ACCMM 
CLl 


1 OKjtVV 

ISGW** 

CL2 

CL3 


=F 0.00065 
± 0.00098 


=F 0.00111 
=F 0.00116 


Total : Lepton Samp 


e 




± 0.0015 


± 0.0035 


r ^ 


1.548 ps 


± 0.032 


± 0.00005 


± 0.00004 


TB± 


1.653 ps 


± 0.028 


=F 0.00002 


=F 0.00003 


TB° 


1.493 ps 


± 0.062 


=F 0.00025 


± 0.00001 


^-^baryon 


1.208 ps 


± 0.051 


± 0.00004 


=F 0.00002 


1 \ 
\~^hadron / 

tb±/t b o 


1.564 ps 


± 0.014 


± 0.00005 


± 0.00000 


1.062 


± 0.029 


± 0.00016 


± 0.00001 


a 

fb— bar yon 


0.115 


± 0.020 


=F 0.00006 


± 0.00010 


/bo 


0.117 


± 0.030 


± 0.00051 


± 0.00001 




0.472 ps' 1 


± 0.017 


=F 0.00003 


± 0.00008 




0.1177 


± 0.0055 


± 0.00082 


± 0.00001 


Total : Mixing 






± 0.0010 


± 0.0001 


Misidentified e 




see text 


± 0.00011 


± 0.00021 


Converted gammas in e sample 


± 10 % 


=F 0.00019 


=F 0.00050 


Misidentified 




see text 


± 0.00025 


± 0.00085 


/i from n, K decay 




± 15% 


± 0.00027 


± 0.00092 


background asymmetry 


± 40 % 


=F 0.00076 


± 0.00421 


Pt reweight of background 


see text 


=F 0.00007 


± 0.00021 


Energy flow correction 


see text 


=F 0.00009 


=F 0.00020 



Total : Lepton identification and px measurement ± 0.0009 ± 0.0044 



b- tag tuning 

Merging for N int = 2 

Jet charge stat 

Jet charge BKG subtraction 


see text 
see text 
see text 
see text 


=F 0.00009 
+ 0.00061 
=F 0.00052 
=F 0.00069 


± 0.00028 
- 0.00082 
± 0.00070 
± 0.00084 


Total : b- tag and jet charge calibra 


tion 


± 0.0011 


± 0.0014 


Total 




± 0.0023 


± 0.0058 



Table 7: Different systematics in the x 2 fit of the 1993-1995 DELPHI lepton sample at 
= 91.22 GeV . The systematic errors at = 89.43 GeV were estimated to be ± 
0.0024 for and± O.OO46 for A^ B , and at ^/i = 92.99 GeV the corresponding values 
were ± 0.0024 and ± 0.0070. 
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central values and variation were taken from [13] and [16] . The b and c decay model and 
their associated changes were taken from [15]. The variation considered for the lepton 
identification was described in Section 3. 

Mixing 

The combined effects of the uncertainties quoted in Table 7 for the parameters having 
a direct effect on the B°B° mixing description give a precision of ± 0.014 on f B o and ± 
0.005 on x [13] following the method described in [14]. To take into account correctly 
the impact of the different sources of uncertainty on the B°B° mixing description, each 
of these measurements was varied within its error. 

Background Asymmetry 

It should be noted that, while the observed "Background Asymmetry" systematic in 
Ap B comes from uds events, for , the main source is the charge correlation between 
a fake lepton and the initial quark in b events themselves. Even if the charge correlation 
between the fake leptons and the initial quark of the corresponding event has been taken 
from the simulation, this correlation can be studied in the data using the Qe x Q opp 
observable. Such charge correlation in b events is visible for example in the left plot of 
Figure 2. In this case the ~ 10 % of the difference in the amount of background between 
the two extreme Qi x Q opp bins, originates from b — > c — > s — > K~, where the K ± is 
misidentified as a muon and behaves like a right sign lepton. 

From these studies a conservative change of ±20% in the charge correlation, corre- 
sponding to ±40% variation in the background asymmetry, has been considered for the 
systematics. Such a change increases by 1.8 the % 2 of the data/simulation comparison 
computed in the left plot of Figure 2 corresponding to 6-tag events enriched in fake lep- 
tons by a.p,pr cut. The same comparison done in the anti-6-tag bin (I), enriched in fake 
leptons from uds events, gives an increase of 1.3 in the \ 2 ■ 

The Ap c B systematics are dominated by the contribution of the background asymmetry. 
This underlines the difficulty to separate the c events from the other flavours (see Fig- 
ure 3). The size of this systematic follows the variation of the asymmetries with y/s and 
explains most of the 50% increase in the A^ B systematics between y/s = 89.43 GeV and 
y/s = 92.99 GeV . 

Energy flow, p T reconstruction 

Due to a slightly worse energy reconstruction in the data, a 1-2 % shift in the jet energy 
distribution between data and simulation has been observed. This difference could have 
different effects on the asymmetries depending on its exact source ( overall correction 
or sub-sample of charged/neutral track correction ). The different possible sources were 
considered and the biggest effect observed was taken as the systematic error. 
In the anti-6-tagged sample the shape of the pt distribution of the lepton candidate was 
not correctly described by the simulation. This effect is known to be common to all 
tracks from hadronisation in the tuned DELPHI simulation [8]. The full size of the cor- 
rection estimated in the anti-6-tagged sample, was considered as a systematic error. It 
corresponds to changes ~ ±5% of the number of misidentified leptons as a function of pt- 

6-tagging 

To take into account the effect of changes in the fraction r q , defined in Equation (5), the 
6-tagging corrections, 5q\ were recomputed for each of the changes quoted in Table 7. 
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For this reason all the quoted systematics also include possible variations induced by 
changes in the 6-tagging tuning. 

The systematic named b-tag tuning in Table 7, corresponds to the effects of the finite 
simulation statistics used to estimate the sample composition in the different 6-tagged 
intervals and to the sensitivity to the correlation pf 1 as described in Section 4. 

The full difference between the results obtained for the two considered N int choices 
(see Section 4) is quoted as Merging for N int = 2. 

Jet Charge 

A jet charge tuning was performed for each computation of the systematics. For this 
reason, possible systematic errors in the jet charge tuning arising from the variation of a 
given parameter defined in Table 7 are included in the systematic errors of this parameter. 

The systematic named Jet charge stat in Table 7, corresponds to the effect of the finite 
statistics used to estimate o~Q. etl , f and Qj et i- The systematic named Jet charge BKG 
subtraction in Table 7, is associated to the uncertainty on the jet charge description of 
non b events (see Section 5). 

8 QCD Corrections to the measured asymmetries 

The QCD corrections applied to the asymmetries were obtained following the prescrip- 
tion given in [17]. This approach takes into account changes in these corrections due to 
experimental bias, like the suppression of events with an energetic gluon induced by the 
cut on the momentum of the selected leptons. The simulation sample, with an enlarged 
asymmetry 12 to improve the statistical precision of the study, was used to estimate such 
a bias. The relative change in the corrections due to experimental bias was estimated for 
this analysis to be 0.58 ± 0.08 for and 0.42 ± 0.12 for A^ B . These scale factors were 
applied to the theoretical QCD corrections 13 and give the following QCD corrections : 
^noQCM = A xy(! _ Cx) with Cb = 0>0 205 ± 0.0046 and C c = 0.0172 ± 0.0057 . 

9 Conclusion 

The heavy flavour asymmetry measurements presented in this paper, obtained with the 
1993-1995 DELPHI data, can be combined with the 1991-1992 DELPHI measurements 
of Ape and Ap B using leptons [3]. All asymmetry measurements were QCD corrected 
before averaging and the 1991-1992 DELPHI measurements were corrected to the same 
inputs (branching ratios and mixing) as the ones used in this paper. 

Taking into account the correlations between the different systematic sources, the 
combined results are: 



Abb 
^FB 


= 0.067 


± 


0.022 


(stat) 


± 


0.002 


(syst) 


at 


V5 = 


89.43 GeV 


Abb 
^FB 


= 0.1004 


± 


0.0056 


(stat) 


± 


0.0025 


(syst) 


at 




91.26 GeV 


Abb 
^FB 


= 0.112 


± 


0.018 


(stat) 


± 


0.002 


(syst) 


at 


v^ = 


92.99 GeV 


ACC 

^FB 


= 0.031 


± 


0.035 


(stat) 


± 


0.005 


(syst) 


at 


v^ = 


89.43 GeV 


A cc 
^FB 


= 0.0631 


± 


0.0093 


(stat) 


± 


0.0065 


(syst) 


at 


v^ = 


91.26 GeV 


A cc 
^FB 


= 0.110 


± 


0.028 


(stat) 


± 


0.007 


(syst) 


at 


v^ = 


92.99 GeV 



12 A value of 0.73 was used, slightly smaller than the maximal asymmetry allowed (0.75) to avoid the boundary problem 
and consequent asymmetric errors in the result of the fit 

i3 c had,T = 0354 ± o 0063 and c had,T = q Q413 ± 0.0063 as recommended in [16]. 
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[GeV] 

Figure 5: Measured values of and A^ B at different y/s compared to the Standard 
Model predictions for m t = 175 GeV/c 2 and mu = 300 GeV/c 2 . (To avoid the overlap 
of the results at peak-2 and peak+2, the and A™ B points have been shifted respectively 
by +0.02 GeV and -0.02 GeV .) 

Figure 5 shows the energy dependence of A^ and Ap B compared to the Standard Model 
prediction. 

Following the general procedure described in [15], these results have been corrected 
to the Z pole; the energy shift from y/s to rriz (for y/s = 89.43 GeV : +0.0391, 
+0.0997; for </b = 91.26 GeV : -0.0013,-0.0034 ; for y/s = 92.99 GeV : -0.0260, 
—0.0664), the effects of the initial state radiation (+0.0041, +0.0104) and 7 exchange and 
7/ Z interference (—0.0003,-0.0008) have been corrected by adding the quoted numbers 
respectively to A^ and A^ B . The averages of the pole asymmetries obtained after these 
corrections are : 

A ^ = 0.1021 ± 0.0052 (stat) ± 0.0024 (syst) 
A% = 0.0728 ± 0.0086 (stat) ± 0.0063 (syst) 

in agreement with other DELPHI and LEP measurements [2,18,19]. 

The total correlation between Ap B and A^ B is +7% , with a correlation of +22% and 
-36% for the statistical and systematic errors respectively. 

The effective value of the weak mixing angle derived from these measurements is 

sin 2 9^ s = 0.23170 ± 0.00097. 
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